The dynamics of a qubit coupled strongly with a quantum oscillator 
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The dynamics of a qubit coupled with a quantum oscillator is re-studied in the region of strong 
coupling. The non-degenerate perturbation is added to the usual degenerate one and new results 
are given. 
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I. INTRODUCTION 

A qubit coupled with a quantum oscillator is the sim- 
plest model in quantum system, yet is an abundance of 
applications both in theoretical and practical aspects. Its 
dynamics is described by the Rabi Hamiltonian [l| . How- 
ever, its analytical solutions is still unavailable. Many 
approximations are applied. Among them, the rotation 
wave approximation (RWA) [2| and adiabatic approxima- 
tion (AA) are major ones. Under weak coupling and near 
resonant condition, it has been solved analytically by the 
method of RWA by discarding the non-energy-conserving 
terms, and is also called the Jaynes-Cumming model in 
quantum optics0- In the opposite condition, that is, the 
strong coupling and large detuning, it is investigate by 
the adiabatic approximation (AA) and generalized RWA 
3-0. The A A method assumed large detuning and first 
only treat the quantum oscillator influenced by the qubit 
through the coupling, then solve the whole system by 
treating the qubit as a perturbation. Further assumption 
of the qubit must be made for its solutions in AA method, 
that is why sometimes called it adiabatic approximation. 
The adiabatic assumption requires that the qubit as a 
perturbation does not have influence on the transition of 
the quantum oscillator from its different number states 
and discards these transition terms concerning the dis- 
placed number states and solved it accordingly. 

Because recent achievements in circuit QED make the 
strong or even the ultrastrong coup ling regime of light 
matter interaction achievable j3l| - [1 C| . this in turn renews 
the theoretical interesting in qubit under adiabatic ap- 
proximation. Further extensions are made by replacing 
the qubit by two qubit and N-qubits or linear coupling 
by non-hnear ones, etc 

In the paper, we will restudy the qubit coupled strongly 
with quantum oscillator and focus on the rationality of 
the adiabatic assumption. Though physically the low en- 
ergy of the qubit can not promotes the transition of the 
oscillator to its different displaced number states, these 
omitting terms sometimes are larger than the keeped 
terms of the transition of the same number states dis- 
placed differently. So we will make a scrutiny into the 



adiabatic approximation and a generalized approxima- 
tion is given and its new results are discussed. Section 2 
reviews the adiabatic approximation in the study of the 
qubit and quantum oscillator coupled system and section 
3 gives new approximation. Final section will discuss and 
conclude the results. 



II. THE BRIEF INTRODUCTION OF THE 
ADIABATIC APPROXIMATION 

The Rabi Hamiltonian for a qubit interacting with a 
harmonic oscillator is Q , Q , 

H — —fkjjoaz + hjoa) a + h(3[a + a^)axi (1) 

where ctz are the usual Pauli matrices with ~ 
a+ + a-. 

In the paper, we do not consider the RWA case where 
the non-rotation-terms a~'"cr+, aa- are neglected. Under 
large detuning and strong coupling condition, the qubit 
is regarded as a perturbation. The so-called free hamil- 
tonian is 



whose solutions are 



(2) 



i7o|^,n,m) = [N - l3^)nw\N^,m), (3) 
|7V„,m) = |m)|iV,„), m = +, (4) 

where |+), |— ) are the eigenvectors of the operator 
with ax\+) = \+), cr^|-) = -|->, and \Nm) are the 
displaced number states for the quantum oscillator, that 
is 



\N±)^D{tP)\N), D{/3) 



^ g/3(at-a)^ 



(5) 
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|+), |— ) are not | t)j \ 1) ^ which are the eigenvectors of 
the operator cr^. 

The Hamiltonian describes the quantum oscillator in- 
fluenced by force from its coupling with with qubit, and 
its solutions show that the force of the qubit acting on 
the oscillator has make the oscillator displace its equilib- 
rium position according to the states of the qubit, just 
as shown in the displaced number states Itpn.m) in its 
solutions. 
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The displaced number states \N±) are not completely 
orthogonal due to the fact 



(iVzp|iV±) 7^ with m ^ m' . 



(6) 



This will result in the complex situation when including 
the perturbation term ^HuJodz of the qubit to be treated. 
This term ^hwQCTz of the qubit will produces the transi- 
tion of the oscillator from its various eigenstates \N±). 
From the theoretical view, the large detuning (wq ^ uj) 
will make it reasonable to omit its transition among dif- 
ferent numbers states. This means that the qubit mainly 
results in transition of the oscillator from the same num- 
ber states which are displaced differently induced by the 
corresponding states of the qubit , that is, the transi- 
tions of the form |A^±) — )• |iV=p). The other transitions 
{\N±) ^ \N'^),N ^ N' ,) will be omitted from the con- 
sideration. This approximation is called adiabatic. In 
this way, the whole Hamiltonian becomes block-diagonal 
as 



/i?o 

Hi 



H = huj 



••• \ 









Hn 



(7) 



with the 2x2 matrices defined as 
„ _ f N-I3^ Oat 



(8) 



where 
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cxp(-2/32)L^(4/32). 



(9) 



L jv (x) in the above equation are the Laguerre polynomi- 
als. The new solutions are easy to obtain: 



4=(|iV+,+>±|iV-,-». (11) 



All eigen- values Epf _j_ arc influenced by the parameters 
/3, A^, so are the eigenvectors _(_). 

The eigenvectors _|_) are complete and orthogonal 
basis for the composite system of the qubit and the quan- 
tum oscillator. The evolution of the qubit will definitely 
depends on the initial states of the composite system. 
Here we only give some examples for showing the dynam- 
ics of the qubit. When the initial states of the composite 
system are i?°(0) = |\E'Ar.+), the possibility of the qubit 
remaining in its initial state are 



P(t) = COS^{nN *Ujt). 



(12) 



The general initial state will evolute by the principle of 
superposition, sec details in Ref.jl]. 




FIG. 1: The dynamics of the system depends on the param- 
eters n, j3. The thick, solid and dashed lines correspond to 
the figures of P{t) with A = 1, /3 = 0.2, iV = 4, /3 = 0.2, 
N = 1, /3 = 0.7 respectively. 



III. THE PERTURBATION RESTUDIED IN 
DETAILS 



As stated before, the adiabatic approximation omits 
various terms of the form 

|A±) ^ |A^^), N^N' 

under the physical ground that uiq <^ uj. Though physi- 
cally the qubit seems not able to promote the oscillator 
to go from one displaced number state to a new one with 
different number as its energy gap \n — n'\hujQ is far more 
less that hw of the transition |A^+) |A^), N ^ N' 
or |A''_) \N'_), N ^ N' , however, this is not true 
for the transition of the type |A+) — > |Ai), n ^ n' 
or |A^_) |A|), N ^ N' . Detailed quantities study 
also support the conclusion, and shows that the per- 
turbation terms {N±\^fkjJoaz\Nj.) depend non-linearly 
on the coupling strength /? and the parameters N,N'. 
Given an arbitrarily strong coupling /?, it is generally 
true that {N^\^hu}oCTz\N±) will be larger than most 
(10) of {N'\^hLOoaz\N±), N ^ N' . However, there are 
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still several numbers N = Ni, i = 1,2,--- , A''^ that 
{N^\^hajoaz\N±) are much larger than {N^\^Hujo(Tz\N±) 
. Quantity study shows the number A^^ is finite and 
depends both the coupling strength /3 and the num- 
ber A^, so it will change accordingly. For example, 
Fig.® shows (13+, -|-|ifiajocrz|10-, — ) is greater than 
(10+, -|-|^/iwoCTz|10_, — ) when /3 = 0.7. So the adiabatic 
approximation must add these terms to be more applica- 
ble. 

There are two kinds of perturbations in quantum the- 
ory. The perturbation to the degenerate system and that 
to the non-degenerate ones. Let's reconsider the pertur- 
bation HluqCz to the free Hamiltonian Hq. The eigen- 
states are |A^+, |A^-,-), A^ = 0, 1,2,---. The states 
with the same numbers n are degenerate, while other 
states with different numbers are not. So the perturba- 
tion ^hujQaz will relate the two kinds ones in quantum 
theory. Obviously, the adiabatic perturbation is the one 



FIG. 2: The comparison of the various perturbation terms 
{10-\^hjJoaz\N\.) . The sohd and dashed hues correspond to 
the figures for (10_ | i&Joa^ |iV4) with /3 = 0.2, 0.7 respec- 
tively. The independent variable is A'''. 



FIG. 3: The comparison of the various perturbation terms 
{0-\N!^). The solid , thick and dashed lines of the figures 
represent for (0_|Af^} with /? = 0.2, 0.5, 0.7 respectively. 
The independent variable is A'^'. 



only applied to the degenerate cases of the free Hamilto- 
nian. So, it is not complete as there is the non-degenerate 
perturbation to need to be added. 

Some observation is given here for convenience: when- 
ever there are perturbations of combination of both the 
degenerate one and non-degenerate one, one must first 
treat the degenerate one, then treats the non-degenerate 
one. The reverse order can not work, as the condi- 
tion of perturbation is small relative the original en- 
ergy gap of the free Hamiltonian can not be met due 
to its degeneration. Further more, for the second per- 
turbation, there are two bases to use: one is from 
the eigen-states of the the free hamiltonian, the other 
comes from the first degenerate perturbation. In our 
case, they are |iV+,-h), |iV_,-), N = 0,1,2,--- , and 
= 0, 1, 2, • • • , i = +, -. Of course, the basis 
\E^^ = 0, 1, 2, • • • , i = — are more favorable 

for the calculation of the second perturbation. 

By the use of the basis iV = 0, 1, 2, • • • , the 

dynamic equation for the composite systems is 

H\En,+)=En,+ \En,+), H\En,-) ^ EN,-\EN,-).{li) 

Then the corresponding eigen- vectors and eigen-energies 
are 



\E 



\E, 



K+) +E«/I^"+) (14) 
1=1 

\E%.^) + EK/I^" +> + bNj\El^)), (15) 



7=1 



E 



N,+ 



E 



N.- 



= huj 



^ 4(7V-/) ^2 



(1,6) 



Nk 



(TV - /?2 - n^) + 



(1,7) 



where 



^N,I 



'^N,I 



^N,I ^ 



(l + (-l)^-0(/+|jV_) c.o 
4(iV - /) cj ' 

(l-(-l)~-^)(/+|iV_)c.o 



A{N - I) 



(18) 
(19) 



The summation includes those terms of (7V^ \^1tjjJq(Jz\N±) 
that are comparable with the term fijv = 
{Nzi^\^hujQ(j z\N±) . See Appendix for details. 

IV. DISCUSSION AND CONCLUSION 

From comparison of the results ([T5|) .(|T6 |) -(fT7 | with the 
old ones ([TT|) . pUj) . we could see the justification of the 
adiabatic approximation. The parameter ^ is small, the 
new eigenstates ([T5|) differ from the old one (fTT]) only on 
first order of Though the new eigenvalues only has 
some extra terms in the second order of — , the old one 
(|10|) is already up to the first order of Some terms 
{N'j^\^hjjjQaz\N±) may be comparable or even larger than 
{N^\^hijOQ(Jz\N±) ^ discarding these terms to make the 
whole Hamiltonian block-diagonal generally is feasible 
under the large detuning condition. However, because 
the term f2„ depends non-linearly on the parameter /3, in 
case of it being zero by some chosen /3, that is at the crit- 
ical points, we consider it is reasonable to refer the adi- 
abatic approximation as the new results ([T5|) . ([T5|) -([T7 1) . 

Here we give an example. Consider the composite sys- 
tem is initially in state *(0) = |-h)|0+) = |0+,+). The 
quantum oscillator is in its displaced vacuum state or 
coherent state |0+) = £)(— /3)|0). Suppose the coupling 
strength /? = 0.2, we will have 



1 



'|0+,- 



El^ = hu{-l3^ + no), El_ = fi^(-/32 - no), (21) 

UJO_ 

2lo 



e-'^'\0-,-)], (20) 

'-%), (21) 
(0-|0+> (22) 
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as the state of the system at time t by the old adia- 
batic approximation method. From Fig. ([3]), we obtain 
that (0_|0+) = 0.98099, (0_|1+) = 0.19604, (0_|2+) = 
0.0277242, (0_|3+) = 0.00320132 and all others much 
smaller than 0.00320132. So, in the new results of 
Eqs. (fT5| . we have 



\E'> 



-) + (421^^2,+) + ^o+il^^i" -)) + b+,\El_)l23) 



If we further choose — = 0.3, then 

I,! ' 



\E. 



0, + / 



Similarly, we have 



\E, 



0,- 



K-) 



0.0029406|£;5'_> 



0.00294061^^ 



(24) 



(25) 



The interesting things of the above equations are that 
|£'o,+) combines the two states \Eq _^_), and |i?o.-) 

combines the two states \Eq _), \Ei _^_) with opposite co- 
efficient. Of course, the coefficient in front of \Ei _) or 
\Ei ^) very small. The change of the eigenvalues is ig- 
norablc in this case, that is, Epf^± sa Ef^ ^. Therefore, 



= e '~ 



'|0+, +) becomes 



0.0029406 
0.0029406 



(e- 



(26) 
(.27) 



some possibility to evolve into its displaced one photon 
states. The result may be of some use for the application 
of qubits coupled with a quantum oscillator. 

In summary, concerning the system of a qubit coupled 
strongly to a quantum oscillator, the adiabatic approxi- 
mation method is restudied and is extended or modified 
to the second order of ^ for the eigenvalues and the 
first order for the eigenstates. Though the modification 
is small and the old adiabatic approximation is justified, 
the new results will have some application. We could 
also extend the study to the system of iV— qubits cou- 
pled with a quantum oscillator, where the spectrum of 



the free Hamiltonian Hn is El, 



0, ±1, 
AE° = 



N 



0,1, 

.2«2 



N.i 



The smallest energy gap is 



huj{l~Am f3 ). For example. Am 



3, the condition that the qubits energy Hujq is much more 
smaller than the transition energy gap AE = fia;(l — 3/3^) 
can not be met for strong coupling (large (3). So it is 
plausible the inclusion of the iV-qubits will not induce 
the the transition of different displaced number states, 
that is, the adiabatic approximation might be put into 
doubt in this case. Our study in the paper must be used 
to see that the adiabatic approximation can still be jus- 
tified. This will be investigated in ref.[ll|. 
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APPENDIX 



Eg. pT]) has the extra term, its last two parts, compared 
tith the old one ([22]) . So, even the quantum oscillator 
is initially in its displaced vacuum states, it will have 



By the use of the perturbation method in quantum 
mechanics, we could solve Eg.p^ by the following 



= +5^6+ = \E%,_)+Y.''-^j\E':+)+T.^Nj\El-), (28) 

1=1 I—l 1=1 i=l 



where 



'■N.I 



^0 /j^O 






{El- 


^^\E%,+) 


_ Wo 




{N 


-l)huj 


A{N - 


I)UJ 


{EU 


^^\E%J 


Wo 




(N 


- i)huj 


4(iV- 


I)UJ 


(^?.+ 


^^\E%,+) 


Wo 






- i)huj 


4(iV- 


I)UJ 



(-(/+|7V_) + (/_|7V+)), (31) 



^Nj = ',t : - ,,,, \^^ (-{I+\N.) (/_|7V+)), (32) 
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and 



^ Nk 1 

En,+ = + -hL^oY^'^NjUN+li-) + {N-\i+)) + -hLJoY.^h(-iN+\i-) + 



1=1 



1=1 



K 



|(7V+|/_)-(7V_|/+)|^ 



1 1 



1=1 



1=1 



, (33) 



■|(Ar^|J_)-(iV_|/+)p + ^ 



16(iV - 



H(7V+|J-) + (iV-|J+)|t3d) 



In all the above equations, the summations all the above results concerning the coefficients a^j, etc, 
mean the sum of the terms involved the quantities 



{N'^\^fujjQaz\N±) that can not be ignorable. Usually, 
there are just several such terms need to calculate. 
Because of the fact Q 



(M+|7V_)=e-2^-(-/3) 




-M-N 



could be simplified as Eqs. ((T5|) . Please note that L^{x) 
in above is the associated Laguerrc function. 

The simple relation among the coefficients 
j,ajf J J, bjf J shows that the all upper lim- 
its of summation in Eqs (fT5 |) - (fT7| could be put to the 
(4/3^), M > iV(35) same as Nk which depends on both N, I. So, the results 



and {M+\N^) being real, there are the relations 
{M+\N^) = (-l)^-^^(M_|iV+) 

and 

(A/_|7V+)=(-l)^-^^(7V_|M+), 



could be concisely written as Eqs (IT5|) - (IT71) . 
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